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ABSTRACT 

We model the inner knot of the Crab Nebula as a synchrotron emission coming from the 
non-spherical MHD termination shock of relativistic pulsar wind. The post-shock flow 
is mildly relativistic; as a result the Doppler-beaming has a strong impact on the shock 
appearance. The model can reproduce the knot location, size, elongation, brightness 
distribution, luminosity and polarization provided the effective magnetization of the 
section of the pulsar wind producing the knot is low, a < 1. In the striped wind model, 
this implies that the striped zone is rather wide, with the magnetic inclination angle 
of the Crab pulsar > 45°; this agrees with the previous model-dependent estimate 
based on the gamma-ray emission of the pulsar. We conclude that the tiny knot is 
indeed a bright spot on the surface of a quasi-stationary magnetic relativistic shock 
and that this shock is a site of efficient particle acceleration. On the other hand, the 
deduced low magnetization of the knot plasma implies that this is an unlikely site 
for the Crab’s gamma-ray flares, if they are related to the fast relativistic magnetic 
reconnection events. 

Key words: ISM: supernova remnants - MHD - shock waves - gamma-rays: theory 
- radiation mechanisms: non-thermal - relativity - pulsars: individual: Crab 


1 INTRODUCTION 

The Crab pulsar and its pulsar wind nebula (PWN) re¬ 
main prime targets for high energy astrophysical research. 
In many ways, the current models of Active Galactic Nuclei 
and Gamma Ray Bursts are based on what we have learned 
from the studies of the Crab. The recent detection of flares 


from the Crab Nebula by AGILE and Fermi satellites (Ta- 
ri et al.||2011 Abdo et al.|2011 ) have brought this object 


into the “focal point” once again. Their extreme properties 
seem impossible to explain within the standard theories of 
non-thermal particle acceleration and require their overhaul 
with important implications to high energy astrophysics 
in ge neral (e.g. |Lyutikov||2010| Clausen-Brown & Lyutikov 


2012[|Cerutti, Uzdensky &; Begelman|2012||Lyubarsky|2012[ 

Buhler fc Blandford|2014| 

In the MHD models of the Crab Nebula, the super-fast- 
magnetosonic relativistic wind of the Crab pulsar terminates 
at a reverse shock ( Rees &; Gunn|[l974 Kennel & Coroniti 


1984). However, finding the shock in the images of the Crab 


Nebula has not been a straight-forward matter - there seem 
to be no sharp feature which can be undoubtedly identified 
with the shock surface. In their seminal paper, |Kennel 
Coroniti (1984) discuss the under-luminous region hosting 


the Crab pulsar and surrounded by the optical wisps as an 
indicator of the shock presence. After the discovery of the 
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inner X-ray ring by Chandra ( Weisskopf et al.||2000 Hester 
et al.||2002 ), the ring is often referred to as the termination 
shock and yet this feature looks much more like a collec¬ 
tion of knots than a smooth surface. A new twist in the 
story has come with the recent PIC simulations which show 
the shock particle acceleration is highly inefficient in even 


relatively weakly magnetized relativistic plasma (Sironi & 
Spitkovsky 2009, 2011b). These results make one doubt that 


the shock can be visible at all. On the other hand, the wind 
from an oblique rotator should have the so-called striped 
zone where the orientation of magnetic field alternated on 
the scale of the pulsar period. The magnetic energy associ¬ 
ated with these stripes can be dissipated at the termination 
shock and converted into the energy of the wind particles 
(Lyubarsky 2003 Sironi & Spitkovsky 2011a|a Amano & 


Kirk|2013 Petri &; Lyubarsky||2007 ). 


Given the highly anisotropic nature of the wind, the 
termination shock is squashed along the polar direction and 
can be highly oblique with respect to the upstream flow 
(Ly ubarsky| [2002). Downstream of the shock, the flow can 
still be relativistic and its emission subject to strong Doppler 
beaming. The computer simulations of the Crab nebula and 
its radiation (Komissarov & Lyu barsky| [2004) revealed the 
presence of a very bright compact feature in the synthetic 
synchrotron maps, highly reminiscent of the HST knot 1 of 
the Crab Nebula located very close to the pulsar (also called 
the inner knot, Hester et al. 1995]). (In these simulations, the 
termination shock was treated as source of synchrotron elec¬ 
trons with power-law energy spectrum, which then were car¬ 
ried out into the nebula by the shocked wind plasma.) This 
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2 Lyutikov et al. 


et al. 2009) and 3D (Porth, Komissarov & Keppens 2014) 


feature, confirmed in the later more advanced 2D (Camus where g\ — Bi /p[ is the magnetization parameter of the 

wind, B' and p are the comoving values of the magnetic field 
and the rest-mass density of plasma respectively. This is a 
monotonic function increasing from x(0) = 1/3 to x(+°o) = 
1. For (j\ 1, one can use the approximation 


simulations, is associated with the location at the termina¬ 
tion shock where the shocked plasma flows in the direction 
of the fiducial observer and thus strongly Doppler-boosted. 
Komissarov & Lyutikov (2011) argued that given the short 


synchrotron life-time of the high energy electrons compared 
to the dynamical time-scale of the shock, the knot can be 
the main source of the gamma-ray emission from the Nebula 
at 10-100 MeV. 

Recently, a targeted multi-wavelength study of the 


Crab’s inner knot has been conducted by Rudy et al. (2015) 
in order to check if it shows any activity correlated with the 
gamma-ray flares. Although no such correlation has been 
found, the optical data reveal the structure and temporal 
evolution of the knot with unprecedented detail. In this pa¬ 
per, we investigate if the data are consistent with the MHD- 
shock model of the knot using simple analytical and semi- 
analytical tools. In particular, we combine the theoretical 
shape of the shock with the oblique shock jumps in order 
to obtain the Doppler-beaming of the post shock emission 
and use this to determine the location, the shape and the 
brightness distribution of the knot. 


2 GEOMETRY OF THE TERMINATION 
SHOCK 

At the location of the termination shock, the magnetic field 
of the pulsar wind has the form of loops centered on the pul¬ 
sar’s rotational axis. The wind’s termination shock is also 
symmetric with respect to the rotational axis and hence the 
magnetic field is parallel to the shock surface. The pulsar 
wind is not spherical - its luminosity per unit solid angle 
increases with the polar angle measured from the pulsar 
rotational axis. As a result, the termination shock is not 
spherical and the radial stream lines of the wind are gener¬ 
ally not normal to the shock surface - locally the shock is 
oblique. In addition, the pulsar wind is ultra-relativistic and 
its thermal pressure is negligibly small. The corresponding 
shock equations have been analyzed in |Komissarov &; Lyu-| 
tikov ( 2Q11|; Lyutikov, Balsara h Matthews (2012); see also 
Appendix [A] Here we summarize their results using the no¬ 
tation introduced in Komissarov & Lyutikov (2011). 


We differentiate the flow parameters upstream and 
downstream of the shock using indices “1” and “2” respec¬ 
tively. Denote as 5 the angle between the velocity vector 
and the shock surface (the angle of attack). Then in the 
observer’s frame 

tan52 = xt an ^i (1) 

and for the Lorentz factor of the flow 

r 2 = Ti [l + r? sin 2 5! (1 - X 2 )] ~ 1/2 , (2) 

where x = v n 2 /v n i is the ratio of the normal velocity com¬ 
ponents. For a strong shock, T 2 <C Ti and the last equation 
reduces to 


r 2 = (l-x 2 r 1/2 csc,5i. 


(3) 


Assuming <5i 1/ri and using the ratio of specific heats 
7 = 4/3, [Komissarov &; Lyu tikov] ( |20111 ) obtained 


X = 


1 + 2(71 + y/l6(7i + 16(71 + 1 
6(1 + ( 71 ) 


X ~ 1 2(71 ' 

Using Eq. § we find that for = 0 

r 2 = —^ 7 = esc 5i 

2v^ 

and for ai 1 
— ^/aiCScSi . 

The deflection AS = Si — S 2 is given by 

A _ tan £i(l — x) 

tan AS = - -^. 

1 + xtan Si 

It reaches the maximum value of 

tan(A5 max ) — ~ — at tan 5i — x ^ 

2 y/x 

For g\ — 0 this gives AS max = 7t/ 6 at £1 = 7 r/ 3 , whereas for 
ai 1 one has A Smax — 1/4c7i at S± = 7 r/ 4 . 

_eF 

The total pressure p 2 = p 2 H— ^ downstream of the 
shock is 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


P 2 = (1 - x)(F/c) sin 2 <5i, 


( 10 ) 


where F is the upstream total energy flux density along the 
flow velocity (see Appendix |A|). For ai = 0, this yields 

P 2 = ~— sin Si (11) 

o c 

whereas for ai 1, 

1 F . 2r / 10 x 

P 2 — - -sin Si . (12) 

2(7 1 c 

One can see that for the same energy flux the post-shock 
pressure is significantly reduced compared to the purely hy¬ 
dro case. 

Since the shock is driven into the wind by the pressure 
inside the nebula, p n , which is approximately uniform in the 
nebula due to its slow expansion, we replace p 2 with constant 
p n , which makes our approach similar to the Kompaneets 
approximation Kom paneets|1 960). This approximation was 
already used by Lyubarsky (2002), to determine the shape of 


the termination shock for a weakly magnetized wind. It less 
clear if the approximation can hold well for the polar sec¬ 
tion of the shock where the magnetization and the Lorentz 
factor of the postshock flow can be very high. This makes 
terms other than the total pressure potentially important 
in the transverse force balance. This is already seen in the 
numerical simulations with moderate wind magnetization, 
where the magnetic hoop stress leads to compression of the 


polar region (Porth, Komissarov & Keppens 2014). More¬ 


over, these simulations show that the polar flow is highly 
variable. Keeping these in mind, we shell still shell proceed 
exploring the models based on the assumptiopn p 2 = p n = 
const. 

If the function R(9) gives the spherical radius as a func¬ 
tion of the polar angle on the shock surface then 


(4) 


tan<51 = f w 


(13) 
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Figure 1. Top panel : Shape of the termination shock for a uni¬ 
form distribution of cr\ . The curves correspond to o\ =0,1 and 
10, from top to bottom. Middle panel: Shape of the termination 
shock for the non-uniform distribution of o\ corresponding to a 
striped wind with a = 7t/4 and f(0) = sin 2 6. The lines cor¬ 
respond to (Jo = 0,1,10 and 10 3 , from top to bottom. Bottom 
panel: Shape of the termination shock for the non-uniform dis¬ 
tribution of (7i corresponding to a striped wind cro = 10 3 and 
f(6) = sin n 6. The lines correspond to a = 80°, 60°, 45° and 20°, 
from top to bottom. The solid lines correspond to n = 2 and the 
dashed lines to n = 4. 


For an axisymmetric radial wind, its energy flux can be writ¬ 
ten as F — Lof{0)/4irR 2 , where f(6) describes the wind 
anisotropy. We will consider only / = sin n 6 , where n — 2 


for the monopole model of the pulsar magnetosphere (Bo- 


govalov 1999). Recently, Philippov, Spitkovsky X Cerutti 


(2015) argued for n — 4, based on their numerical simula¬ 
tions of pulsar magnetospheres. Substituting the expressions 
for sinJi and F into Eq. (10), we obtain the shock-shape 
equation 


r, 2 +r, =W} ( i - x) - (i4) 

Finally, we introduce the characteristic length scale of the 
problem Ro = Lq f (0)/4 tt cp n and arrive to the dimensionless 
equation 


X' 2 + X 2 = f(0)(l - x), 


(15) 


where X = R/Rq. (This is the modified version of our orig¬ 
inal Eq.3.) The appropriate boundary condition is 


X(0) = 0. 


(16) 


When the shock terminates the striped part of the pul- 



19 

Figure 2. Magnetization of the striped wind zone after the dis¬ 
sipation of stripes for the magnetic inclination angles a m = 60 
and 45 degrees. 


sar wind, the shock solution is modified due to the dissi¬ 
pation of the magnetic energy associated with the stripes. 

) have shown that the shock solution is ac- 
as that for the unstriped how where the 
energy of stripes is already converted into the bulk kinetic 
energy of the wind particles. Thus, as long as the shock so¬ 
lution is concerned it does not matter where the dissipation 
occurs, in the wind or at the shock. The magnetization of 
the wind that has lost its stripes can be found as 


Lyubarsky (2003 


tually the same 


cri — cro 

Xoc(0) 


(17) 

1 + CT 0 (1 ~Xa(d)) ’ 


where 




Xa{0) = 

j (2M0)K-1) 2 , 

0m < 0 < 7r/2 

0< 0m 

(18) 

and 

cosda = 

tan 0 jji 



(K omissarov|201 3). In these equations, 0m = 7r/2 — a, where 
a is the pulsar’s magnetic inclination angle, is the polar 
angle of the boundary separating the unstriped polar section 
of the pulsar wind from its equatorial striped zone and cro 
is the original magnetization of the striped wind. Figure [2] 
shows the wind magnetization after the dissipation of its 
stripes for cro — 100 and <a m = 60° and 45° degrees. The 
most interesting feature of these solutions is the rapid drop 
of a at the boundary of the striped zone. 

Eq. (15) is integrated numerically. Due to its singularity 
at 0 — 0, we use its asymptotic analytic solution 


x = 2(1 -x(<ro)) e %+i 

77/ T 2 


in order to move away from the origin. 

As a start, we consider the case of uniform cri, where it 
does not depend on the polar angle. This corresponds to the 
case of aligned rotator, a = 0, where cr\ — cro everywhere. 
The top panel of Figure [l] shows the solutions for a i = 0,1 
and 10. As one can see, for higher g\ the shock is located 
closer to the pulsar. This is in agreement with the earlier 
results by Kennel X Coroniti (1984). In fact, the curves differ 
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4 Lyutikov et al. 


only by the scaling factor -y/l — x( <T i)> as follows from Eq. 

d- 

When the variation of <7i due to the existence of the 


striped wind zone (see Eq. (17)) is taken into account, the 
variation of the shock size is less dramatic. In the middle 
panel of Figure [l] we show the solutions for a = 7r/4 and 
f(9) — sin 2 6 , corresponding to different values of the mag¬ 
netization parameter cr 0 . As cr 0 increases, the shock still be¬ 
comes more compact, but as <7o oo the dependence be¬ 
comes very weak and the shock approaches some asymptotic 
shape. Such a turn is clearly connected with the existence 
of the striped wind section where ( 0 ) becomes insensitive 
to (Jo: 

(1-XcW ( } 

The bottom panel of Figure^illustrates the dependence 
of the shock shape on the magnetic inclination angle a for 
do = 10 3 and f(0) — sin n 0 with n — 2,4. As one can see, 
the shock becomes more compact as a decreases. This is 
expected, as for a = 0, the case of uniform magnetization 
with <ji = 10 3 is recovered and in this case the shock size 
rapidly decreases with <to. However, even for a = 20° the 
equatorial radius of the shock is still much larger than that in 
the limiting case of a = 0. For n=2, the total wind power is 
L w = (10/15)Lo, whereas for n=4 it is L w = (8/15)Lo- Since 
it is more interesting to compare the results corresponding 
to the same wind power, we re-scale the n=4 solution of Eq. 
(15) by the factor y/5/2. In the bottom panel of Figure [l] the 
n=2 solutions are shown as solid lines and the n=4 solutions 
as dashed lines. The difference between the two groups is not 
large, particularly for a < 45°. 

Figure [3] zooms into the inner region of the middle panel 
of Figure where the shock exhibits a noticeable break. 
The origin of this break is easy to understand. At 0 < 0 m , 
the magnetization <ri = <jo is constant. Hence, the shock 
curve is a miniature version of that of pure hydro shock 
and 12). At 6 — 6m, &i{6) rapidly drops leading 


11 


(see Eqs 

to higher wind “ram” pressure and the shock shoots out 
almost radially until the “ram” pressure approaches that of 
the nebula. This interpretation suggests that for high cro, 
the shape of the equatorial part of the termination shock is 
independent on that in highly magnetized polar section. 

The low ram pressure of the termination shock in the 
high-cr polar region and the rapid drop of a± around 6m 
suggest that, as far as the equatorial part of the shock is 
concerned, one can ignore the presence of the polar section of 
the wind altogether. In this approximation, the appropriate 


boundary condition for Eq. (15) is 


X(6m) = 0 . 


( 20 ) 


Figure [4] compares this approximate solution with the orig¬ 
inal one for a = 20°, the case in Figure [l] with the largest 
unstriped sector. In this case, the difference between the so¬ 
lution is expected to be most profound. Yet, as one can see in 
this figure, it is still rather small. This result is particularly 
welcome as one expects to see significant deviation from the 
uniform pressure distribution of the shocked plasma in the 
polar region where the high-sigma post-shock flow remains 
supersonic. The exact details of the flow in this region should 
not matter much. 

Strictly speaking, our analysis shows that there is no 
well defined unique shape of the termination shock which 
can be used to predict the emission properperties of Crab’s 
inner knot. On the other hand, the dependence on the wind 



r 

Figure 3. Shape of the termination shock near the origin. The 
solutions correspond to the model with cro = 10 3 and f{6) = 
sin 2 6. As in the middle panel of Figurejl] the magnetic inclination 
angle a = 80°, 60°, 45° and 20°, from top to bottom. 



Figure 4. Shape of the termination shock for ao = 10 3 and 
f(6) = sin 2 0 and a = 20°. The solid line shows the original 
solution and the dashed one the approximate one for the trun¬ 
cated wind and the dash-dotted one the solution corresponding 
tO (71 = cro- 


parameters is not that strong. With the exception of very 
small magnetic inclination angle, the shock shape is approx¬ 
imately the same as found for the weakly-magnetized wind 


by Lyubarsky (2002). For this reason, we will use this shape 


for the rest of our paper. After small additional rescaling, 
the shock shape in this case is described by 


X' 2 + X 2 = f{0). 


( 21 ) 


With f(Q) = sin n 0, the asymptotic solutions of Eq. (21) 


X ~ e — (l - J-6 2 

2 V 48 

for n — 2 and 

x - e ~x(i- 2 

3 V 90 


for n — 4. The corresponding angles of attack are 

Si ~ 0/2 

and 


( 22 ) 

(23) 

(24) 


<5i ~ e/3 (25) 

respectively. Figure [5] illustrates how the termination shock 
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Figure 5. View of the polar region of the termination shock for 
f(6) = sin 2 6 and the viewing angle 0 o b = 60°; the pulsar position 
is shown by the dot. 
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Figure 6. The observed optical emissivity of the post-shock flow 
in the 3D RMHD simulations of |Porth, Komissarov &; Keppensj 
|201 4). The arrowed lines are the instantaneous stream lines. The 
dashed line is the line of view and the blue curves show the regions 
of enhanced observed emissivity. 


with n = 2 appears to a distant observer for the viewing 
angle of 60 degrees to the symmetry axis. 


3 ESTIMATES OF BASIC PARAMETERS. 


As the shocked plasma expands and slows down, its observed 
emissivity drops. Fig ure |6| based the results of 3D RMHD 
simulations by Porth, Komissarov <V Keppens (2014), illus¬ 


trates this behaviour. One can see a relatively thin layer of 
enhanced emissivity just above the shock surface. Its thick¬ 
ness is approximately one third of its distance from the line 
connecting the origin (pulsar) and the observer. The main 
reason for the drop of the emissivity with the dinstance from 
the shock is the reduction of the Doppler beaming. 


Komissarov & Lyutikov (2011) estimated some of the 


the knot parameters in the shock model, assuming that they 
are determined by the Doppler boosting of the emission from 
the shocked plasma. In their calculations, they assumed that 
the velocity of the plasma is parallel to the shock surface. 
Here we do a more careful and extended analysis. 

Assuming a small size of the knot, we first ignore vari¬ 
ations of the proper emissivity across the knot. In this 


case, the observed synchrotron emissivity is (e.g., Lyutikov, 
Pariev & Blandford 2003, see also © ' 


(X x> 2+(p ~ 1)/2 \B\ | (p+1)/2 . 


(26) 


where p is the spectral index of the electron energy spec¬ 
trum, B'j_ is the normal to the line of sight component of 
magnetic field in the fluid frame and V = r -1 (l — i>cos<a) -1 
is the Doppler factor. Even if the magnetic field strength is 
constant over the knot, B' ± may still vary significantly across 
the knot due to the relativistic aberration of light. However, 
along the symmetry axis in the plane of the sky B' ± — B ', 
and it is only the Doppler factor that matters. 


3.1 Low a at the knot location 

Based on Equatiorj26] one can immediately rule out g\ > 
1 for the termination shock at the location of inner knot. 
The key observational data here is the clear separation of 
the knot from the pulsar (Rudy et al. 2015). This shows 
that that the beaming angle ad is smaller compared to the 
deflection angle AS of streamlines at the shock. Defining ad 
as the angle at which D 2 +( p-1 )/ 2 reduces the factor of two, 


we find that for the observed spectral index p « 2.5 

ad ~ 2iy 

Using the maximum value for the deflection angle and T 2 
for high < 7 i (Equations |7| and [8]) we find that 

■^4 ~ 2y / oTsin5i y/2 dq , 

L\0 

where we used S 1 = 7r/4 as the angle of attack with maximal 
deflection. For the case of Si <C 1, we find that 

Both results show that for <ri > 1 one has ad > A S and 
hence the pulsar has to be embedded into the knot, in con¬ 
tradiction with the observations. Figure [ 7 ] shows the ratio of 
angles as a function of <j\ for £1 = 7t/4 and 7r/20. One can 
see, that the dependence of Si is rather week. Using [ 7 ] and 
[8]one can show that for <j\ <C 1 

~ 0 7 
A5 ’ 

and thus the knot size is comparable with the separation 
from the pulsar. This conclusion does not depend on the 
shape of the termination shock and thus very robust. 

Since is expected to be low only in the striped-wind 
zone, this allows us to conclude that the magnetic inclination 
angle a > 90° — 0 o b ~ 30°, where 0 o b ~ 60° the observed 
angle between the line of sight and the rotational axis of the 
Crab pulsar (N g &; Romani| |2004). Based on this result we 
focus in the rest of the paper on the case of low o \. 


3.2 Separation from the pulsar 

The brightness peak of the knot corresponds to the point 
where the deflected streamline points directly towards the 
observer. The polar angle of this point 0k — 0 o b — A S (see 
Figure®. Using Eq. © in the limit of small angles and 
Eq. (|24|) , we find that 

A$«(l-x)tfi«(l-x)| (27) 

and 

dk = ■ (28) 

The angular separation between this point and the pulsar in 
the plane of the sky is 
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Figure 8. Geometry of the shock in the plane defined by the location of an observer, the pulsar and a point on the shock surface. 



lo Sio CT i 

Figure 7. The ratio of the beaming and deflection angles. 


ipk 



(29) 


where D is the distance to the pulsar. Denoting as Rt s and 
ipts the linear equatorial radius of the termination shock and 
its angular size in the plane of the sky respectively, 


i>k 

'Ipts 

= it AS - 

JXts 

(30) 


For cr <C 1, one has x~l/3 and hence 


AS 

~ -0 o b , Ok = -0 o b • 

(31) 


Now Rk = 0.28 and 'ipts = HAip p ~ 7.4". This is approxi¬ 
mately equal to the radius of the Crab’s halo ( |Hester et al. 


1995) and almost twice as small compared to the radius of its 
X-ray ring. For the shock shape function f(6) = sin 4 6 , one 


obtains Rk = R s ( 0k ) — 0.14i2o, Rts = R s ( 7t/2) 0.70i?o 
and ipts = ZO'ipp ~ 13". Thus, given the uncertainty of the 
shock shape, the theory and observations are quite consis¬ 
tent in the limit of low <j\. 


3.3 Transverse size 


The full half-brightness transverse size of the knot can be 
estimated as 


=2 a h ^. 


(32) 


where an is the angle between the line of sight and the veloc¬ 
ity vector at the point on the shock, with the same position 
on the shock-defining curve as the center of the knot, where 
the emissivity is reduced by the factor of two because of the 
Doppler effect and the relativistic aberration of light. Thus, 
we have 


A^_l _ 2 an 

'ipk AS 


(33) 


The observed synchrotron emissivity is given by Eq. 
(26). Provided the knot size is small, one can assume that the 
magnetic field is uniform and write B' ± = B' cos a , where 
the a is the angle between the stream line and the line of 
sight in the fluid frame. The relativistic aberration of light 
gives 


cos a — v 
1 — v cos a 


Substituting this into Eq. (26), we find that 


e u oc (1 — ^cosa) ^ +2 ^(cosa — v )^ p+1 ^ 2 . 


(34) 


(35) 


Approximating v ~ 1 — 1/2T| and cos a ~ 1 — <a 2 /2, this 
reads 

e u oc (1 + x 2 )~^ p+2 \l — x 2 )^ p+1 ^ 2 , (36) 
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cr=0 <T= 1 cr=10 



Figure 9. Synthetic emission maps of the knot for cr% = 0, 1, 10 (from left to right) and f(6) = sin 2 6. In all cases, the distances are in 
arcseconds, the pulsar is located at the origin and the emission peak is at 0.7 7/ from the pulsar. The peak emissivity is normalized to 
unity and the third contour corresponds to one half of the peak value. For high cr > 1 the inner knot is both very broad and elongated 
in radial direction, which is inconsistent with the observations. The contours start from 0.9 of the peak value and decrease by the factor 
of thereafter. 



Figure 10. Sensitivity of the synthetic images to the shock model. Left panel: shock shape f(6) = sin 2 6 and uniform proper emissivity; 
Center panel: shock shape f(0) = sin 4 0 and uniform proper emissivity; Right panel: shock shape f(0) = sin 2 0 and proper emissivity 
scaling with the spherical radius as oc R ~ 2 . In all the cases a\ = 0. The contours start from 0.9 of the peak value and decrease by the 
factor of y/2 thereafter. 


where x — T 2 <a. For the observed p = 2.6, this equals to 
one half for x « 0.33 and. Thus ah ~ 0.33/T2 and Eq. (33) 
reads 


A'if± ^ 0.66 
ifk T2A(5 


(37) 


Substituting into the last equation the expressions © and 
© in the approximation of small Si, we finally obtain 


A ^l 


0.66 


1 + X 
i - X 


1/2 


(38) 


This is a monotonically increasing function of x an d has 
the absolute minimum value (A'ip±/'i/jk)min — 0.9 reached 
for G\ — 0 (x = 1/3). For a i = 10 this gives A^>±/'ifk = 
4.3. Observational measurements of the knot parameters are 
complicated by its small size and proximity to the bright 
Crab pulsar. Depending of the method used, the transverse 
size of the knot in HST images varies from Ax/j± to 

0.56 7/ , whereas 'ifk ~ 0.65" (Ru dy et al.|2015 ). This rules out 
high a i and favors ai < 1 once more. The synthetic images 
of the knot presented in Section [4] give a somewhat smaller 
size compared to what follows from Eq. (38). 


4 SYNTHETIC IMAGES 

In this section we construct two-dimensional “images of the 
knot”. Obviously, in order to obtain the brightness distri¬ 
bution we need to integrate the emissivity along the line of 
sight. However, from the shock geometry we can only con¬ 
clude how it is distributed over the shock surface. Therefore, 
in this section we start by constructing images of this surface 
and later study the effects of finite thickness of the emitting 
layer. We expect a longer geometrical length of the emitting 
region along the line tangent to the shock surface. This fac¬ 
tor would make the knot more compact along the symmetry 
axis in the image. On the other hand, the finite thickness 
of the emitting layer would tend to increase the knot size in 
this direction. 

For all images presented in the paper we use 0 o b = 7r/3. 


4.1 Emissivity maps 

Given the shock shape and the “upstream” magnetization 
parameter cj\ one can determine the post-shock flow direc¬ 
tion and its Lorentz factor, as well as the angle \i — 7t/2 — a 
between the line of sight and the magnetic field in the fluid 
frame. These allow us to compute the purely geometrical 
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Figure 12. Emissivity maps for the striped-wind shock shape. Top row: a = 7r/4. Bottom row: a = 7r/3. Left column: cro = 10, n = 2; 
Middle column: cro = 100, n = 2; Right column: cro = 10, n = 4. The contours start from 0.9 of the peak value and decrease by the factor 
of y/2 thereafter. 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 


Figure 11. Shock contour for obliqueness of 7r/3, 7 t/4, 7t/6 (top 
to bottom) for wind magnetization of cro = 10 and f(0) = sin 2 0. 
The hydro shock solution is indicated dashed black. The shock- 
contours are colored with the Doppler factor cubed. Red contour 
describes the los to the Pulsar for the assumed inclination angle 
$ob — 60°. Right-hand panel is a zoom-in. One can see the double¬ 
humped structure of the termination shock, leading to a emission 
closely aligned with the los in the case of a = 7r/4. 


component of the synchrotron emissivity over the shock sur¬ 
face. Namely, Eq|26| gives us that 

€u oc Cgeom = V 2+(<p ~ 1)/2 1 s in fi '| (p+1)/2 . (39) 

Next, we project this distribution of e ge0 m on the plane of 
the sky. The main contribution to the knot emission comes 
from the closest to the observer section of the shock surface. 
Due to the non-spherical shock geometry, the line of sight 
may intersect this section twice. In this case, we sum the 
contributions from both these points. Next we rescale the 
image so that the maximum is located at 0.7 arcsec from 
the pulsar. 

As an illustration, Figure]!)] shows the results for the 
shock shape described by Eq[21| with n — 2 and constant 
magnetization parameter <j\ — 0,1,10. One can see that 
only in the case g\ — 0 thereis a clear separation between 
the knot and the pulsar, in full agreement with the results 
of Sec [3] The plots also confirm the conclusion of Sec|3] that 
high <71 models result in radially elongated elongated images, 
which is in conflict with the observations. 

In Figure |10| we compare the results for n — 2 and 
n — 4. One can see that the difference between the model 
is not dramatic - in the n — 4 model, the knot is a little 
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Figure 13. Knot images resulting from emitting region of finite thickness (top) and corresponding emissivity distributions (bottom). 
The dashed line shows the line-of-sight and the white dotted line traces the assumed shape of the termination shock. The contours start 
from 0.9 of the peak brightness and decrease by the factor of \/2 thereafter. 


bit more tangentially elongated. The proper shock emissiv¬ 
ity may reduce with the distance from the pulsar reflecting 
the reduced wind power. To probe the inportance of this fac¬ 
tor, we also considered the model where the shock emissivity 
scales as R~ 2 e ge0 m ~ the results are shown in the right panel 
of Figure [lO] One can see that the knot becomes significantly 
more compact and less elongated in the radial direction. This 
image is closer to those of the Crab’s inner knot, which is 
approximately 2:1 in size (tangential over radial), while its 
separation from the pulsar, ~ .65 arcsec is much larger than 
it’s radial width (0.15 arcsec in the HST image and some¬ 
what larger, 0.35 arcsec in the Keck image), ([Rudy et al. 
20l5 ). 

Finally, we have also explored the case of the shock 
shape of striped wind described by Eq. [14] with <j\ varying 
according to Eq[l7| In all models described here, the shock 
emissivity is e = R~ 2 e ge0 m■ 

For large magnetic inclination angle a > 0 o b , the shock 
shape is very similar to the our “standard” one. Moreover, 
(Ti <C 1 and hence the images are not much different from 
those shown in Figure |10| For small magnetic inclination 
angle a < 0 o b the knot emission comes from the inner lobe 
of the shock, where the shock shape is exactly the same as 
the standard one. However, a i is very high now, leading to 
images which are in stark conflict with the observations (like 
the one in the left panel of Fig. [9|. 

The most interesting is the case with a ~ 0 o b , where 
the knot emission comes from the transitional section of 
the shock where a± varies rapidly and the shock surface is 
closely aligned with the line of sight (see Figure [IT]). In Fig- 


12 we show the results for a = 7r/3 and 7t/4. One can 


see that for a = 7r/3, a the images can show a very high de¬ 
gree of elongation in the transverse direction. However, for 
a = 7 t / 4 , this elongation is no longer seen. The results con¬ 
firm our expectation that the differences between standard 
shock shape (Eq[l5f and that of the striped wind ( |21| are 
not dramatic and the conclusions based on the models with 
standard shape are quite robust. 


Rudy et al. (2015) also pointed out that the Crab’s 


knot could be a bit convex away from the pulsar (the “smily 
face”). In the synthetic synchrotron maps presented in this 
section, the distant side of the knot has a sharp edge slightly 
convex the other way. This feature reflects the curvature of 
the folding edge of the shock surface projection onto the 
plane of the sky. However, in the images based on the RMHD 
numerical simulations of PWN the edge curvature is washed 
away (Porth, Komissarov & Kep pens||2014| . In these simu¬ 
lations, the emission comes from a layer of finite thickness 
downstream of the termination shock (see Figure [6|. As we 
show next, this can be an important factor in determining 
the detailed shape of the knot. 


4.2 Brightness distribution 

In order to probe the effect of finite thickness of the emitting 
layer, we need a model of volume emissivity away from the 
shock surface. Our starting point is the emissivity on the 
shock surface, which we assume to be 

e o {R,0,<j>) oc R~ 2 V 2+{p ~ 1)/2 \sin f j,'\ {p+1)/2 . (40) 
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Figure 14. Polarization of the shock integral emission in the 
case of zero thickness as a function of the particle spectral in¬ 
dex p. Solid curves correspond to the integration over the whole 
shock surface for models with a i = 1, 0.1,0 (top to bottom) and 
constant rest-frame emissivity. The dashed line shows the result 
for spherical surface obtained in fLyutikov, Pariev &; Blandford| 
(2003}. 


The emissivity outside of the surface is then modeled as 

e(R,0) = eo(R,0,<l>)6e(R,e), (41) 

where the function S e (R : 0) provides spreading about the 
surface. We choose it to be 

S e (R,e) = exp(-\R-Rs(e)\/(e r Rs(0))) (42) 

for R > R s (9) and 

5 e (R,0)tcmh(-4\R- R s (0)\/(e r R s (0))) + 1 (43) 

for R < R s (0), where R = R s (0) is the shock radius. The 
parameter e r controles the relative thickness parameter. The 
factor of 4 in the argument of tanh provides much faster drop 
of emissivity in the direction towards the pulsar. 

In Figure [T3] illustrates the results obtained for the 
shock shape parameter n = 2 with <j\ — 0. As one can 
see, the “frown” turns into a “smile” already for e r = 0.1. 
At this point, the emissivity is still a sharp layer attached to 
the shock. Notwithstanding the ad-hoc nature of this sim¬ 
ple model, the results indicate that the shape of the knot 
is very sensitive to the downstream flow. We hence suggest 
that a modeling of the knot’s shape must take into account 
the flow in the post-shock layer. 


5 POLARIZATION 


Given the velocity field and the assumed magnetic structure 
of the flow (azimuthal) we can also calculate the polarization 
signature. In order to do this properly, the relativistic aber¬ 
ration of light has to be taken into account Lyutikov, Pariev] 
&; Blandford (2003). In our case, the calculations are slightly 
different due to the different geometry of the problem. The 
details can be found in Appendix |B] 

We start with the case where emission come only from 
the shock surface. Figure [14] shows the degree of polarization 
for the total flux coming from the shock with n — 2 and 
proper emissivity scaling as R~ 2 . (the results for n — 40 
look very similar.) One can see that for 0 < a± < 1 the 
degree of polarization varies only slightly. For the observed 
value of p = 2.6, we obtain II « 50%. For magnetization 
<ti > 1, the polarization signal nearly coincides with that 


Figure 15. Polarization fraction in the model with finite width 
of the emitting layer ( g\ = 0 ). The width is scaled with the local 
shock radius: 10% - dashed, 20% - dotted and 50% solid. Even 
for a unrealistically wide emitting region of 50%, the observed 
polarization degree of ~ 60% still cannot be reproduced. 


for a spherical shock (Lyutikov, Pariev & Blandford 2003) 
- in this case the flow deflection angle A 5 is small and its 
speed is highly relativistic. 

Thus, our model predicts high polarization of the knot 
emission in agreement the observations, but the predicted 
value is still somewhat lower than the observed one of II ~ 
60% ( Moran et al.|2013 ). In order to understand the reason 
we carried out additional polarization calculations. 

To check if the finite thickness of the emitting layer 
can effect the polarization, we carried out calculations with 
the same volume emissivity model as in Sec |4.2| The results 
are presented in Figure p~5] which shows that the degree of 
polarization remains largely unchanged - the changes are 
quite small - when the thickness is increased from 10% to 
50% of the local shock radius, the polarization increases by 
merely 2%. This leads us to conclude that a finite extent 
of the emitting region alone is unlikely to explain the high 
observed degree of polarization. 

In the observations of Moran et al. (2013), the knot 
polarization was measured in a very localized area with an 
aperture radius of 0.15 // . However according to the data by 


Rudy et al. (2015), the transverse FWHM size of the knot 
is ~ 0.32 and FWRMS size is « 0.56" and hence the aper¬ 
ture used in Moran et al. (2013) captures only the bright 
inner part of the knot. This can be significant because the 
depolarization of total flux in our calculations is caused by 
the gradual rotation of the polarization vector across the 
knot, which is illustrated in Figure [16] Thus, a smaller area 
of integration would give a higher polarization degree. 

In order to investigate this effect we carried out addi¬ 
tional calculations where the integration over the azimuthal 
angle was limited to the interval (—fa,+ <j>b). To determine 
a reasonable range for fa, we recall that the knot emissiv¬ 
ity decreases by the factor of two from its peak value for 
streamlines making angle ah ~ 0.33/F2 ~ 0.14 to the one 
leading to the peak (see Section pL3} . 

In Fig. [IT] we show the results of integration for fa = 
0.1, 0.25 and 0.5. One can see that for all these value the 
polarization is significantly higher compared to what we ob¬ 
tained previously. In fact, for fa — 0.1 the flux polariza¬ 
tion degree almost coincides with the theoretical maximum 
in uniform magnetic field. Deviations from the exact axial 
symmetry of our model will naturally reduce the polariza¬ 
tion degree. 
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Figure 16. Image of the knot (cri = 0) and its polarization vectors (E-field) for f(0) = sin 2 0, (left) and f(6) = sin 4 6, (right). At 
each location the polarization degree corresponds to the theoretical maximum for the synchrotron emission. The rotation of polarization 
vectors across the image results in depolarization of the integral emission. The contours start from 0.9 of the peak value and decrease by 
the factor of y/2 thereafter. 



Figure IT. Polarization signal coming from the area of the shock 
limited by \<f\ < 0.1, 0.25, 0.5. (At the peak of emission the ab¬ 
scissa of 0.15 arcseconds corresponds to </> « 0.1.) The dashed 
line is the theoretical maximum for the synchrotron emission in 
uniform magnetic field. 


6 OTHER PROPERTIES 
6.1 Energetics 

Let us estimate the energy flux intercepted by the region 
producing the wind and compare it with the observed lu¬ 
minosity. To calculate the solid angle occupied by the knot, 
we recall that for ai < 1 the part of the flow which con¬ 
tributes to the knot emission occupies AO « 0.3 with a 
A« 0.66/T2 ~ 0.3 in the case ai C 1 (see Sec[3]). In 
the same limit, the deflection angle AS « ffi(l — x) ~ 0k/ 3, 
where 0k is the coordinate of knot center. This gives 0k — 
0 o b — AS ~ 36U/4 ~ 7r/ 4 and Si « 7t/8. The wind lumi¬ 
nosity per unit solid angle L(0) = L w sin 2 0 , where L w is 
the total wind power. Thus, the energy passing through the 
knot is 


Lknot 


4 x 10~ 3 L U 


= 2 x 10 36 erg/s, 


(44) 


where L w = 5 x 10 38 erg/s is the current spin-down power of 
the Crab pulsar. Given the observed spectrum of the nebula, 
most of this energy is carried out by the electrons emitting 
in the optical band. 

The observations give the isotro pic optical-IR luminos- 
ity of the knot of 1.3 x 10 33 ergs -1 (Rud y et al.|2015 ). Given 
the Doppler beaming angle of I/T 2 with Y 2 ~ — 2.5 


1.6 x 10 31 erg/s, 
5 


the actual total knot luminosity is L 0 b 
implying the radiative efficiency of f ra d ~ 10 

The synchrotron life-time of optical electrons is 

a 2 x 10 5 Sl| /2 days, 


t 


(45) 


where B- 3 is the magnetic field in mGauss. Taking the knot 
size along streamlines of ti c ~ 10 light days as a reasonable 
estimate, the knot crossing time in the fluid frame will be 
ti c ~ 4 days. Hence, the theoretical radiative efficiency of 
the knot f ra d = ti c /t syn « 2 x lO -5 # 3 ^ 2 , which is consistent 
with the observations. 

For comparison, the total isotropic luminosity of bright 
wisps within 10" from the pul sar is about ten ti me that 
of the inner knot, ~ 10 34 erg/s (Hester et al. 19950 Their 
proper motion indicates velocities v ~ 0.6c. Hence the beam¬ 
ing angle in the 0 direction is A0 ^0.6 and in the <f> direc¬ 
tion A (j) « 0.5, slightly smaller due to the anisotropy of the 
synchrotron emissivity in uniform magnetic field. The cor¬ 
responding solid angle is about unity and hence the actual 
wisps luminosity is L 0 b ~ (1/An) x 10 34 erg/s. The luminosity 
emitted in all directions will be higher on average by n/ A<f>, 
yielding L W i sps « 5 x 10 33 erg/s. According to the MHD the¬ 
ory, the wisps are arc-like structures of enhanced magnetic 
field advected in the equatorial direction (Ca mus et al.|2 009 


Porth, Komissarov &; Keppens|2014 ). It is in this direction 


where most of the pulsar wind power is transferred. Hence 
the radiation efficiency of the wisp region 


frad — 


10 “ 


(46) 


which is similar to what we found for the kno10 


6.2 Variability 


Moran et al. (2013); Rudy et al. (2015) discuss the variability 


of the position, size and the luminosity of the inner knot - the 
position fluctuates relative to the mean by approximately 


1 We used the wisp length of 3 ", as stated in ( Hester et al.|l 995t 
for the “thin wisp”, to calculate the ratio of the isotropic wisps 
luminosity to that of the inner knot. 

2 The maximal spectral power in the Crab Nebula comes out 
in UV-soft X-rays, where the radiative times scales of leptons is 
roughly comparable to the age of the nebula. Most of this emission 
comes from the old volume-filling population of particles, and not 
from the freshly accelerated ones close to the termination shock. 
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10% on the time scales of month(s). At the same time the 
overall size of knot correlates with the distance from the 
pulsar, while the luminosity anticorrelates with it. 


The numerical simulations by Camus et al. (2009); 


Porth, Komissarov &; Keppens (2014) show that inner re¬ 
gion of PWN is highly dynamic and the shock surface is 
constantly changing as the result. When the external pres¬ 
sure drops the shock expands and when the pressure in¬ 
creases the shock recedes. The emission of wisps is one ef¬ 
fect of this variability observed in the synthetic synchrotron 
images obtained in the simulations. The other one is the un¬ 
steady behavior of the inner knot, whose position and bright¬ 
ness change in time. In fact, |Porth, Komissarov Sz Keppens| 
(2014) reported an anti-correlation of their synthetic knot 


luminosity with its projected distance from the pulsar. 

In order to understand these results, let us consider 
the simplest model of the shock variability, where the shock 
shape is preserved but its length scale Ro fluctuates. In this 
case, the downstream emissivity is the same function of 6 
and 0 up to a factor depending on Ro . This means that the 
ratio of the knot size to its separation from the pulsar re¬ 
mains unchanged (see Sec[3]), which is in good agreement 
with the HST observations, which give oc '0°‘ 8±013 . 

Regarding the total flux from the knot, we note the 
emissivity e v oc n' 2 B 2 ( ' p+1 ^ 2 , where n the number density 
of emitting particles. The total flux of the knot F v oc e v A, 
where A is the knot area. Since, n 2 oc Rq 2 , A oc Ro, B' 2 oc 
Rq 1 and ipk oc Ro we obtain 

F v oc 'ipl where s = — (p ± l)/2 ; 


(47) 

the same result as stated in Rudy et al. (2015). Thus, the 
shock model is consistent with the observed anti-correlation. 


For the observed spectral index p — 2.6, Eq. (47) gives 


s = 1.8, whereas the HST data suggest a somewhat larger 
value s = 2.39±0.37 (Ru dy et al.|2 015). In reality, the shock 
variability may not be shape-preserving, in which case the 
variablity of its observed emissivity will be more compli¬ 
cated. 

The results of computer simulations show that the shock 
variability is more complicated, with the shock shape chang¬ 
ing as well in response to the external perturbations on the 
scale below Ro (Porth, Komissarov &; Kepp ens| 2014). Thus, 
the predictions based on the model of uniform scaling should 
be considered as rather approximate. 

Since in the MHD theory both the wisp production and 
the knot variability are related to the variations of the shock 
geometry, one would expect approximately the same time- 
scale for both these phenomena^] Although the available ob¬ 
servational data do not cover a sufficiently long period of 


time, they indicate that this may be the case (Hester et al. 
2002 Moran et al.||2013 Rudy et ah]|2015 ). 


6.3 Connection to Crab’s gamma-ray flares 


The discovery of flares from the Crab Nebula ( Tavani et 
aT7||2011 Abdo et al.||2011 ) challenges our understanding of 


particle acceleration in PWNe and possibly in other high 
energy astrophysical sources. 

The short life-time of gamma-ray emitting electrons 
means that if they are accelerated at the termination shock 
then the gamma-ray emitting region is a thin layer above the 


Knot variability may occur on shorter scales ( jLyntikov, Balsara| 
Sz Matt hews 1 20 12| . 


shock where the flow Lorentz factor is still high and hence its 
emission is subject to the Doppler beaming. Komissar ov V| 
Lyutikov (2011) used this to argue that most of the observed 


gamma-ray emission of the Crab Nebula may come from the 


inner knot. They and Lyutikov, Balsara V Matthews (2012) 
also speculated that the gamma-ray flares of the Crab Neb¬ 
ula may come from the knot as well and proposed to look 
for correlations between the knot’s optical emission and the 
gamma-ray emission from the nebula. The relativistic post¬ 
shock flow may help to explain the peak frequency of flares 
exceeding the radiation reaction limit of ~ 100 MeV (Lyu- 
tikov|20l0 Komissarov &; Lyutikov|2011 ). Moreover, a blob 
moving through the knot of ~ 10 light days length would be 
observed for the time smaller by 2r|, which can explain the 
short timescales of the Crab’s flares. 

The observed cutoff of the synchrotron spectrum of the 
Crab Nebula at ~ 100 MeV in the persistent Crab Nebula 
emission and especially during the flares, when the cutoff 
energy approached even higher value of ~ 400 MeV, is in 


conflict with slow stochastic acceleration mechanisms (Lyu- 
tikov|20i0 Clausen-Brown V Lyutikov|2012 ). Alternatively, 


the flares may result from linear particle acceleration during 
explosive relativistic magnetic reconnection ( e.g |Lyutikov 


&; Uzdensky [20031 |Lyubarsky||2005| |Lyutikov||2010| |Cerutti 

et al. | |2014| l. Fast and efficient particle acceleration during 
the reconnection requires highly magnetized plasma, a > 1, 
in the flare producing region (e.g., Lyu barsky|2 012). How¬ 
ever, our results show that the knot plasma cannot be that 
highly magnetized. Moreover, the magnetic field in this re¬ 
gion is still expected to be very regular (after the dissipation 
of the small-scale magnetic stripes), namely azimuthal of the 
same orientation (Porth, Komissarov & K eppens|2014 ), and 
hence lacking current sheets required for the reconnection. 
Finally, the coordinated programs of optical observations did 
not reveal anything unusual about the inner knot emission 


during the gamma-ray flares (Weisskopf et al. 


2013 


Rudy 


et al.|2015 ). Thus, we have to admit that the Crab flares are 


unlikely to originate from its inner knot. 

If flares do not come the termination shock then they 
are certainly not connected to the shock particle accelera¬ 
tion mechanism. An explosive magnetic reconnection seems 
to be the only realistic alternative. A favorable location for 
such reconnections would have high magnetization parame¬ 
ter a = B 2 /w > 1, where w is the relativistic enthalpy, as 
this would ensure the relativistic Alfven speed. In addition, 
its magnetic field should be somewhat disordered so that 
thin current sheets may develop. In pulsar wind nebulae, 
such plasma is expected to exist in the polar region down¬ 
stream of the termination shock, which is fed by the un¬ 
striped section of the pulsar wind ( Lyubarsky|2012 Komis- 


sarov|2013 Porth, Komissarov & Keppens 2014 j". 

The current observations do not rule out yet that at en¬ 
ergies below 100 MeV the synchrotron gamma-ray emission 
between flares is coming from the knot. If so, a slow variabil¬ 
ity of the persistent gamma-ray emission at these energies, 
on the time-scale of wisp production, is expected. Additional 
studies are required to clarify this issue. 


7 COMPARISON WITH OTHER STUDIES 

Almost simultaneously with our manuscript, the results of 
an independent study by | Yuan V Blandford (2015) have be¬ 
come publicly available (Some of their results have been out¬ 
lined already in (Ru dy et al.|2015 ).). We agree in the conclu- 
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sion that the observed knot parameters rule out high mag¬ 
netization of the post-shock plasma. However, they could 
not reach a definitive conclusion on the acceptability of the 
shock model even in the low magnetisation regime, pointing 
to a number of difficulties. The main of them concern the 
transverse size of the knot, its shape and polarization. 

(2015) 


For the transverse size, Yuan & Blandford 


claimed that in the basic shock model it exceeds the dis¬ 
tance to the pulsar at least by the factor of 2.8, in conflict 
with the observations. However, this estimate is based on the 
assumption that the emissivity drops by the factor of two 
at an angle = 1/T2 to the velocity vector. In reality, the 
combination of the Doppler beaming with the anisotropy of 
the proper synchrotron emissivity in uniform magnetic field 
leads to a much smaller angle (see Section 3.3). Curiously, 
their synthetic image in Figure 3b shows a much more com¬ 
pact knot, well in line with our results. 

For the polarization of the integral shock emission, they 
obtained a value which is lower compared to that of the 
inner knot as obtained by Moran et al. (2013). In fact, this 
result agrees with our calculations, when the flux integration 
is carried out over the whole shock surface. However, in the 
observations, the polarization is measured only for the bright 
core of the knot (area with an aperture radius of 0.15”). We 
have demonstrated that smaller integration area leads to 
higher polarization degree, allowing a much better fit. 

Finally, Yuan &; Blandford (2015) pointed out that the 
shock model cannot reproduce the “smily” shape of the knot, 
claimed in Rudy et al. (12015), but yields images more rem- 
inicent of a “frown”. This conclusion is base on the model 
where the emission comes only from the shock surface, which 
also leads to a very sharp brightness drop at the distant (rel¬ 
ative to the puslar) edge of the knot. We have shown that 
in models with finite thickness of the emitting layer, these 
features do not survive and the frown can easily turn into 
a smile or even a pout. In fact, the experimentation with 
different geometries of streamlines in the emitting zone by 


Yuan & Blandford (2015) also show rather strong distortions 


of synthetic images. To address such details more advance 
models, based on computer simulations, are required. 


8 CONCLUSION 

In this paper, we have further explored the model of the 
Crab Nebula inner knot as a Doppler-boosted emission from 
the termination shock of the pulsar wind. This model suc¬ 
cessfully explains a number of its observed properties: 

Location: The knot is located on the same side of the pul¬ 
sar as the Crab jet, along the symmetry axis of the inner 
nebula, and on the opposite side as the brighter section of 
the Crab torus. This is a direct consequence of the termina¬ 
tion shock geometry and the Doppler-boosting. 

Size: The knot size is comparable to its separation from 
the pulsar. This also follows from the shock geometry and 
the Doppler-beaming. The anisotropy of the proper syn¬ 
chrotron emissivity, which vanishes along the magnetic field 
direction in combination with the relativistic aberration of 
light is another significant factor. Only models with low 
magnetization of the post-shock flow, with the effective mag¬ 
netization parameter of the wind <j\ < 1 agree with the 
observations. 

Elongation: The knot is elongated in the direction per¬ 
pendicular to the symmetry axis. This is because the knot 


emission comes from the region where the shock surface is 
almost parallel to the line of sight. 

Polarization: The knot polarization degree is high, and 
the electric vector is aligned with the symmetry axis. This 
come due to the fact that the post-shock magnetic field is 
highly ordered in the vicinity of the termination shock and 
azimuthal. In the model, the relativistic aberration of light 
leads to a noticeable rotation of the polarization vector along 
the knot and this prediction could be tested in future polar¬ 
ization observations. Accordingly, the polarization degree of 
the integral knot emission depends on the integration area 
- the bigger the area the smaller the degree is. 

Luminosity: Taking into account Doppler beaming, the 
observed radiative efficiency of the inner knot is consistent 
with efficient particle acceleration at the termination shock 
and the knot’s magnetic field of one milli-Gauss strength, 
which is a reasonable value for the inner Crab Nebula. 

Variability: The knot flux is anti-correlated with its sep¬ 
aration from the pulsar. In the numerical simulations, the 
termination shock is found to be highly unsteady, changing 
its size and shape. As the shock moves away from the pulsar, 
so does the knot region, which leads to lower magnetic field 
and hence lower emissivity. Another outcome of the shock 
variability in the MHD simulations is the emission of wisps 
and hence one expects both the processes to occur on the 
same time-scale, which is consistent with the observations. 


In many cases, the agreement with the observed proper¬ 
ties of the Crab’s inner knot falls short of a perfect fit. Given 
the uncertanties in the shape of the termination shock, 
proper emissity of the shocked plasma and the post-shock 
flow which are present in the model it would be naive to 
expect more. Further investigations of the models, involving 
advanced numerical simulations, are needed to achive this. 

Our results may have a number of important implica¬ 
tions to the astrophysics of relativistic plasma in general and 
that of PWN in particular. They show that the termination 
shock of the relativistic wind from the Crab pulsar is a re¬ 
ality and that this shock is a location of efficient particle 
acceleration. The strong Doppler-beaming of the emission 
from the shock explains why this shock has been so elusive. 
Only the emission from a small patch on the shock sur¬ 
face, the inner knot, is strongly Doppler-boosted and hence 
prominent. For most of the shock, its emission is beamed 
away from the Earth and hence difficult to observe. 

The shock model of the inner knot allows us to constrain 
the parameters of the wind from the Crab pulsar. Taken di¬ 
rectly, the model requires the wind to be particle-dominated, 
<7i < 1 , at least at the polar latitudes of 40° — 60°. How¬ 
ever, in the case of a striped wind, its termination shock 
can mimic that of a low a flow even when the actual wind 
magnetization is extremely high (Lyubarsky; 2003). In this 
context, the magnetic inclination angle of the Crab pulsar 
should be above 45°, which means that most of the Poynting 
flux of the Crab wind is converted into particles, if not in the 
wind itself then at its termination shock (Komissarov] 2013). 
This is in agreement with the results of numerical simula¬ 
tions, which can reproduce the observed properties of the 
inner Crab Nebula extremely well in models with moderate 
wind magnetization (Porth, Komissarov &; Keppens| |2014). 
However, the polar region of a pulsar wind is free of stripes 
and can still inject highly magnetized plasma into its PWN. 

The fact that during the gamma-ray flares of the Crab 
nebula the inner knot does not show any noticeable activ¬ 
ity suggests that the flares occur somewhere else. This is 
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consistent with the fact that any stochastic acceleration 
mechanism is too slow to compete with radiative losses 
and deliver electrons capable of emitting synchrotron pho¬ 
tons of 100 MeV energy. Our conclusion that the inner-knot 
plasma is not highly magnetized also disfavors the knot as a 
site of explosive relativistic magnetic reconnection. To pro¬ 
ceed really fast, the magnetic reconnection has to occur in 

JCL 


magnetically-dominated plasma (Lyutikov 


2010 


Brown &; Lyutikov||2012] JCerutti et al.||2014[ |Lyubarsky 
2012| |Sironi &; Spitkovsky |2014| ). The inner polar region of 


the Crab nebula is the only location where such conditions 
can be met. 
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APPENDIX A: OBLIQUE RELATIVISTIC MHD SHOCKS 


In the shock frame, the fluxes of energy, momentum, rest mass, and magnetic field are continuous across the shock 


(w + B 2 )^ 2 fo = const, (Al) 

2 

(w + B 2 )^ 2 fofo + p + — = const, (A2) 

(w + B 2 )^ 2 fofo = const, (A3) 

PI fix = const, (A4) 

B^yfo — const, (A5) 


where p is the rest mass density, p is the gas pressure, w = pc 2 + kP is the relativistic enthalpy, k — T/(T — 1), where T is 
the adiabatic index, B is the magnetic field as measured in the fluid frame, [3 = v/c , and 7 is the Lorentz factor. We select 
the frame where the velocity vector is in the xy-plane, the magnetic field is parallel to the z-direction, and the shock front 
is parallel to the yz-plane. In what follows we will use subscripts 1 and 2 to denote the upstream an the downstream states 
respectively. 

We assume that the upstream plasma is cold, p\ — 0, and ultrarelativistic, 71 1, that the shock is strong and the 

downstream ratio of specific heats is V 2 = 4/3. Hence fo x ~ sin^i, there #1 is the angle between the velocity vector and the 
shock plane. Denote the wind energy flux in the radial direction as F. Then 


(wi + B\fo\fo x = ( F/c ) sin£1 , 

and Equations |A1 1 and | A2| read 

(w 2 + B 2 )^\fo x = (F/c) sin Ji, (A6) 

(w 2 + B\)^\fo x fo x + P? = ( F / c ) sin2 <*i> (A7) 

where p 2 = P 2 + is the total pressure (Note that we ignore the contribution of the magnetic pressure to the upstream 
momentum flux). Combining the two one finds 

p 2 = (F/c) (sin 2 Si - p x2 sin £1). (A8) 

For a strong shock, 


fox = xfox = X s in ^1 , 


where 

1 + 2cr 1 + 

X 6(1 + <ji) 

( Komissarov &; Lyutikov||2011 ). Hence, 

P‘2 — (1 x)( F / c ) s in 2 Si . 

X is a monotonically decreasing function of ai. For a 1 = 0, one has x — 1/3 and 

~ 2F • 

P 2 = - — sin Si 
3 c 


which is the same as derived in Porth, Komissarov & Keppens (2014). For a 1 1, one has x — 1 — 1/2 <ji and 

1 F . 2 

p 2 - - -sin di . 

2cti c 


(A9) 

(A10) 

(All) 


(A12) 


APPENDIX B: EMISSIVITY CALCULATIONS. 


Let us introduce Cartesian coordinates centered on the pulsar with the z axis aligned with its rotational axis and the 
line of sight parallel to the XOZ plane. In the corresponding basis, the radius vector of a point on the shock surface is 
R s = (sin#cos0,sin0sin0,cos#). The orthogonal projection of this vector into the plane of the sky is 

R S ,_L — Rs (Rs - Hob)rio6 , (-EH) 


where n Q b is a unit vector along the line of sight. In the plane of the sky, we introduce the angular polar coordinates {ipob, (ftob} 
with the origin at the pulsar image and the reference direction given by the orthogonal projection of the rotational axis (see 
Fig.§. Given Eq. ( |Bl| ), the projection of the shock point has the coordinates 


Ipob = ( R s /D ) sin 0 , 

. , R s _l • n y sin </> sin 0 

sm (fiob — | — . 7. 5 

| R s , _l | sinB 

cos 0 = cos 6 cos 0 o b + cos </> sin 0 sin 0 o b 


(B2) 
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and n y = {0,1,0} is the unit vector along the y axis. 

For any proper emissivity, the relativistic Doppler and aberration of light effects ensure that the observed synchrotron 
emissivity 

e„ oc D 2+(p “ 1)/2 |sin y ti , | (p+1)/2 , (B3) 

where p is the particle spectral index and p! = tt/2 — o' is the angle between the magnetic field and the line of sight in the 
fluid frame ( e.g 

At the shock, the post-shock velocity direction is given by the unit vector 


Lyutikov, Pariev & Blandford 2003). 


n f — {cos 4> sin 0 /, sin <j) sin 0 /, cos Of } , 


(B4) 


where Of *= 6 + A 5. Hence the Doppler factor 

D = (r 2 (l-V2(n/-n 0 6 ))r 1 - (B5) 


In the fluid frame, the direction vector of the line of sight is 

n 0 b + r 2 v 2 ^pB2-j-(n 0 (, • V 2 ) — 1^) 
n ob r 2 (1 - (n oi) • v 2 )) 

(Eq. C9 in Lyutikov, Pariev & Blandford 2003) and since the magnetic field is purely azimuthal 
cos// = Uf • n 'ob • 

The unit electric polarization vector (EPV) of synchrotron emission can be found as 
_ n ob X [riff, + n ob x (v 2 x n^)] 

V^l - n ob • v 2 ) 2 - (n 0 • n ofe ) 2 /r| 

(Lyutikov, Pariev & Blandford 2003). The angle \ between this vector and the symmetry axis in the plane of the sky is given 
by 


(B6) 


(B7) 


(B8) 


cosy; = e • (n ob x n y ) , sinx = -(e • n y ). 


(B9) 


According to these equations, it increases in the anti-clockwise direction. Fig. [TB] shows the distribution of the this vector over 
the synthetic image of the knot in the model with f(0) = sin 2 0 and = 0. One can see that the vector is rotating across 
the knot. At each point, the polarization degree is maximal but the polarization of integral emission will be lower due to this 
rotation. 

Due to the mirror symmetry of the image the Stokes parameter U integrates to zero, and the polarization fraction of 
integral emission is 


n = M = p + 1 

I ~ p + 7/3 f t(9,4>)dV 


(B10) 


In models where the emission comes from the shock surface, the emissivity includes the delta-function 6(R — R s (0)), where 
R s (0) is the shock radius. 
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